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Abstract 

We derive the general form of Hawking temperature for Vaidya black hole in the tunneling pictures. This 
kind of black hole is regarded as the description of a more realistic one since it's time dependent decreasing 
mass due to the evaporation process. Clearly, the temperature would be time dependent as our findings. 
We use the semiclassical methods, namely radial null geodesic and complex paths methods. Both methods 
are found to give the same results. Then, we discuss the possible form of corresponding entropy. 
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I. INTRODUCTION 



In 1974, Hawking startled the physics community by proving that black hole evaporates particles 

II a. 

It contradicts with classical general relativistic definition of a black hole, an object that 
nothin g ca n escape from it j^. Hawking's derivation was a quantum field theoretically. With 



Hartle 



22l |. he also derived the black hole temperature by using Feynman path integral which can 



be categorized as semiclassical method, but it still has such mathematical complexity. Semiclassical 
^.Mi.p^„ m were developed in *e pa. decade and «,e. attrac. .an. 



methods for black ho 
attentions, e.g. Refs. 



15|and[l6|. 

There are two ways to perform semi-classical analysis for a black hole radiation. The first is by 

r 1 

the use of radial null geodesic method developed by Parikh and Wilczek [4]. In this method, one 
has to get the expression dr/dt from the radial null geodesic condition, ds"^ = dU = 0, for a metric 
that has the form ds^ = —a (r) dt^ + b{r) dr^ + r'^dQ'^. Then, the obtained expression is used to 
calculate the imaginary part of the action for the process of s-wave emission and relates it to the 
Boltzman factor for emission to get Hawking temperature. 

The second one is called complex paths method was developed by Padmanabhan et.al. 0, Q]- In 
the method, the scalar wave function is determined by the ansatz (j) (r, t) = exp [—iS (r, t)/h] where 
5 (r, t) is the action for a single scalar particle. Inserting this ansatz into the Klein-Gordon equation 
in a gravitational background, one yields an equation for the action S (r, t) which can be solved by 
the Hamilton-Jacobi method. After obtaining the action, one can get the probability for outgoing 
and ingoing particles, P [out] = |(/>out| and P [ in] — |0in| 5 respectively. The Hawking temperature 
can be obtained by using the 'principle of detailed balance' P, Q], -P [out] = exp [— /3i?] P [in] = 
exp [—(3E], since all particles must be absorbed by the black hole. 

In this paper, we consider a more general metric with the mass of black hole is time (t) and radius 
(r) dependent. As generally known, this subject is well described by the Vaidya metric. To get an 
exact {t — r) dependence of the Hawking temperature, we insist to work in using Schwarzschild like 
metric for the Vaidya space time [17[| rather than standard Eddington-Finkelstein metric which is 
widely used for this case. We derive the general form for Hawking temperature both in the radial 
null geodesic and complex paths methods. We arrive only at the general form since one has to get 
the exact (t — r) dependence of mass which of course needs a specific model of the related black 
hole. Until now, even though some models for dynamical black holes had been proposed, e.g. Refs. 



23 and 



24l . but efforts to get better one are still performed. 



The organization of our paper is as follows. In the second section, we will derive the Hawking 
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temperature for a metric with a varying mass by the use of radial null geodesic method. In the 
third section, the Hawking temperature is obtained by the complex paths method, by considering 
only the lowest order expansion of action. At the fourth section, we show the possible expression 
of entropy. In the last section, we give a conclusion for our work. For the rest of this paper, 
we use the unit dimension: Newton constant, light velocity in vacuum, and Boltzman constant, 
G = c = kB = l. 



II. RADIAL NULL GEODESIC METHOD 



but with different 



Previously, a tightly connected work has been performed by T. Clifton 
method of calculations. In this work, we keep that the metric contains the mass of black hole as 
the function of standard radius-time (r — t) coordinate rather than Eddington-Finkelstein radiation 
coordinate {u — v). We start with the metric derived by Farley and D'Eath [l^] for Vaidya space- 
time 

ds^ = - f ^ y (l - ^) de ' dr- + r^dn\ (1) 

\x[m)J \ r J \ r J 

In the above, the black hole mass m varies with time t and radius r, m = m{r,t). The x (m) 

is an arbitrary function of mass and dQ'^ = dO"^ + sin^ Odcp'^ is the metric of unit 2— sphere. We 

will let the arbitrariness of x (m) since the determination of this function depends on the model of 

corresponding black hole. A short discussion of the model will be given at the last section but will 

not the main subject of this paper. The metric ([1]) has the general form 

ds"^ = -F (r, t) dt^ + G (r, t)'^ dr^ + r^dn^, (2) 

where in our case F (r, t) = rn?x (l — 2mr~^) and G {r,t) = 1 — 2mr~^. The general form 

([2]) will simplify our next calculation. Unless for specific purposes, we will write x(m), F{r,t), 
and G (r, t) as x, F, and G respectively for the sake of brevity. 

It turns out that the metric ([T|) has a coordinate singularity at = 2m which of course is time 
dependent. Painleve transformation that is used to remove the coordinate singularity for a metric 
with time-like Killing vector, is also applicable for this analysis. By transforming 



dt^dt-J^^dr (3) 



the metric ^ changes into 



ds"^ = -F (r, t) dt^ + 2F\l ^—^dtdr + dr^ + r^dQ^, (4) 

V FG 



and therefore no coordinate singularity is found. Thus, by such a Painleve transformation, it is 
understandable that, in principle, a coordinate singularity in general relativity can be removed only 
by changing coordinate, without defining new physical condition or theorem. By this thought, it 
would be alright if we still insist in using the metric ([T|) for carrying out a tunneling process. For 
radial null geodesic, ds^ = dil'^ = 0, the differentiation of radius with respect to time can be 
obtained from @ as 

where +(— ) signs denote outgoing (ingoing) radial null geodesies. 

Near the horizon, we can expand the coefficient F and G by the use of Taylor expansion. Since both 
F and G are {t — r) dependent, and we only need their approximation values for short distances 
from a point (horizon), we could apply the Taylor expansion at a fixed time. So, we can write 

Fir,t)\,c^ F'ir,t)\,{r-rh)+0{_ir-rhf)l, (6) 

and 

G (r, t)\, G' (r^, t)\^ (r - + O ((r - r^f) [ , (7) 

By the approximations ([6]) and ([7]) above, the dependence of radius to time in ([5]) can be approached 
by 

dr 1 



dt - 2^'F'irh,t)G>{rH,t){r-TH). (8) 
Now, we discuss the action of outgoing particle through the horizon. In the original work by Parikh 



and Wilczek 



, the imaginary action is written as 

r-out ^^out Vr 



liaS = Im j prdr = Im J J dp',,.dr 



nn '^in 



J J —dr. (9) 

nn 



dr 
dt 



The above expression is due to the Hamilton equation dr/dt = dH/dpr\r where r and pr are 
canonical variables (in this case, the radial component of the radius and the momentum). As a 
reminder, the action of a tunneled particle in a potential barrier higher than the energy of the 



particle itself will be imaginary, pr = y/2m {E — V). Different from discussions of several authors 
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for a static black hole mass, e.g. Refs. I8H151. 1211 and|22, the outgoing particle's energy must be 
time dependent for black holes with varying mass. So, the dH' integration at ^ is for all values 
of outgoing particle's energy, say from zero to +E (t). 

By using approximation ([8]), we can perform the integration ([9|). For dr integration, we can perform 
a contour integration for upper half complex plane to avoid the coordinate singularity r^. The result 
is 

27rE (t) 



ImS" 



Since the tunneling probability is given by F ~ exp 



-|lm5 



(10) 



equalizing it with the Boltzmann 



factor exp [—(3E (t)] for a system with time dependent energy we obtain 



Th 



(11) 



As seen in expression (llip . the Hawking temperature Th is time dependent. We will see it later 
that it is also radius dependent. Partial derivatives of F and G with respect to r are easily to be 
found in the following forms 



F' 



2m 



m' 1 



2m 

r 



mx 



2m 



+ m 



m m 
+ ^ 



and 



G' 



-m + 



m 



Inserting the values of rh = 2m to the last expressions yields Hawking temperature (lllh in a general 
form 

h{l/2-m' irh,t)) 



Th 



(12) 



Attx {rh, t) 

In the above, the m' and x are written has the complete expression as m' {r^, t) and x (r^, t). These 
mean that both m' and x are evaluated at (the radius of event horizon). Our finding in ()12p 
behaves as it is expected for radius independent black hole, where it means that m'{rh,t) = and 
we get Th = Ti/Attx (t) which is close to an expression that had been guessed by Ashtekar long 
time ago, Th = h/Snm (t). 



III. COMPLEX PATHS METHOD 



Massless scalar particles under the gravitational background g^^^ obey the Klein-Gordon equa- 



tion 



--df, [g'^'V^du] </> = 0. 



(13) 



For spherical symmetric black hole, we may reduce our attention only to {t — r) sector in the 
space-time, or in other words, we reduce to two dimensional black hole problems. Later, since 
we are dealing with massless particle described in (jl3p . we can employ the null geodesic condition 
dr/dt = ±V FG for outgoing particle's path. Equation ([13]) under the background metric ([2]) 
simplifies to 

d^<t> + Gf) dt<t^ - \ {F'G + FG') dr4> - fgd^.cP = 0. (14) 

By the standard ansatz for scalar wave function (j){r,t) = exp —^S{r,t) , equation (|14p leads to 
the equation for the action S (r, t) 



Id, 



2l| neglect the terms for n > 1. This 



4<--/-)(^)(f)-/.(i^P)-^(f)>o^ 

Now, our next step is to solve this equation. An approximation method can be applied by expanding 
the action in the order of Planck constant power, 

S (r, t) = So (r, t) + anh^Sn (r, t) , (16) 

for n = 1,2, 3, .... The constant a„ is set to keep all the expansion terms have the action's dimension. 
Taking unit dimensions G = c = ks = 1, otn would have the dimension of [m]^^" which m refers 
to the mass. It is clear that this expansion would lead to a very long equation. Due to the very 
small value of the Planck constant, many authors 
consideration is acceptable, and including higher terms is just adding correction for semi-classical 
derivation of Hawking temperature. By grouping all the terms into the same powers of fi, and take 
the term with lowest order of Planck's constant (zero power), one can write 

FG{drSof -{dtSof = (17) 

Our next task is to find the solution for the last equation which is not too simple since our FG is 
{r — t) dependent. In the standard Hamilton-Jacobi method, 5*0 {r,t) can be written into two parts, 
the time part which has the form of Ft and the radius part 5*0 (r) which is in general a radius 
dependent only. Since our metric coefficients are both radius and time dependent, the standard 

method would not be applicable. We could generalized the method by making an ansatz 

t 

So{r,t) = J F{t')dt' + So{r,t). (18) 



At the first sight, it seems that the ansatz is rather strange, that is S (r, t) and Sq (r, t) are both t 
and r dependent. The term / E {t') dt' is more understandable, since the emitted particle's energy 
is continuum and time dependent. Let see how it works. 
From (jlSp . one can write that 



dtSo{r,t) = E{t)+dtSo{r,t) (19) 

and 

drSo{r,t) = drSo{r,t). (20) 

Since 5*0 (r, t) is t and r dependent, one can write 

^^^^ = drSo (r, t) + dtSo (r, t) ^ . (21) 

Eliminating dt/dr by the use of dr/dt = itV FG^ equation (j2ip can be written as 

"^^^pll. ^ -^dt~S, (r, t) = drSo (r, t) . (22) 
dr \/FG 

Combining the first equation of (jl7p . with equations (jlSp and (j22p where we should note that 
the action equation — (FG) ' dtSo (r, t) = drSo (r, t) is belong to outgoing particle and with the 
radial evolution is dr/dt = VFG, then we could write 

T (FG)-'/' {e {t) + dtSo (r, t)) = ^^^^ T {FG)-'" d^So (r, t) . (23) 

From (f23]) . we can get the exact differentiation of Sq {r,t) 

'-^^ = ^{FG)-'''E{t), (24) 

and the solution of Sq (r, t) can be obtained by integration 

~S,{T,t) = ^E{t) J (25) 

The integration (25) can be evaluated by adopting the value of / (EG) '" dr as in obtaining 
expression (jlOp from Q along with it's approximation method (near horizon Taylor expansion). 
The result for the integration (|25p is 

~ ZTT 

So{r,t) = TE{t)-^. (26) 
The equation gives the complete action 

t 

S{r,t) = j E{t')dt'^E{t)-^. (27) 



The signs + (— ) in expression (j27p refer to the action for ingoing (outgoing) particle. 

Back to our first ansatz for scalar wave function, cj) = exp (r, t) , the wave function for ingoing 

and outgoing massless scalar particle can be read of as 

t 



(r, t) = exp 



(28) 



and 



t>out {r, t) = exp 



-UjE{t')dt'-E{t) 



ITT 



\0 



/TO 



(29) 



respectively. Consequently, from (j28p one can get the ingoing probability of particle as below 

t 



Pin = exp 



^ .Im E (t') dt' + 




(30) 



This ingoing probability must be equal to unity since all particles including the massless one are 
absorbed by the black hole. This consideration gives us the relation 



t 

IiuJe {t') dt' 





ttE (i) 



which leads the outgoing probability 

Pout = exp 



AttE (t) 



(31) 



Finally, to get the Hawking temperature from the outgoingprob ability (jSip . we equate this prob- 
ability expression with exp [—/?£' (t)] which in Refs. |g andlT] is called 'detailed balance' principle. 



It yields 



Th 



h^F' {rh,t) G' {rh,t) 



(32) 



IV. GENERAL FORM OF ENTROPY 



The law of black hole mechanics which expresses the conservation of energy by relating the 
change in black hole mass m to the changes of its entropy angular momentum J, and electric 
charge Q, is given by 2^ 



dm = TndSbh + ^dQ + Q.dJ 



(33) 
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where is the angular velocity and $ is electrostatic potential. By simplifying the black hole with 
no angular velocity and electrical charges, one only has to pay attention on 



dm = TndSbh- (34) 

Inserting our previous result for Th from (jl2p . we an integral form for entropy as 

Stt /■ x{m) 

5,. = -y3-^dm. (35) 

To get an exact value of entropy which is intuitively should be equivalent to dynamical area of 
black hole, we must write x (m), m', and m in their exact forms. This manner would be such an 
interesting work to be pursued further. 



V. SUMMARY 

To summarize our findings above, we can state that we have derived a general expression for 
Hawking temperature (up to the model function x (m) and partial derivative of black hole mass 
with respect to radius) at (fTT]l and (fT2l) in two different semiclassical approaches, null geodesic 
and complex paths methods for a Vaidya black hole. This kind of black hole is regarded as the 
description of a more realistic one since it's time dependent decreasing mass due to the evaporation 
process. Clearly, the temperature would be time dependent as our findings. Both methods are 
found to give the same results. Then, we discuss he possible general form of discussed black hole 
entropy ([35]) . This last finding perhaps can be compared to some models that had been proposed, 



,1 

e.g. Refs. i23 and 



24. 
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